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ABSTRACT 

We show that the recent identification of X-ray emission from SN1998bw is naturally explained as 
synchrotron emission from a shock driven into the wind surrounding the progenitor by a mildly 
relativistic shell ejected by the supernova, the existence of which was inferred earlier from radio 
observations. X-ray observations imply a shell energy E w lO"*^'^ erg, and constrain the initial shell 
velocity /3c and normalized wind mass loss rate, m = (M/10-5MQyr-i)/(ww/10^kms-i), to satisfy 
(3^171 w 10~^'^. The inferred energy is consistent with energy estimates based on radio observations 
provided m « 0.04, in which case radio observations imply (3 « 0.8, consistent with the X-ray 
constraint P^rh w 10~^'^. While X-ray observations allow to determine the parameters characterizing 
the prc-explosion wind and the mildly relativistic shell ejected by SN1998bw, they do not provide 
evidence for existence of an ofF-axis "standard" GRB jet associated with SN1998bw, that may have 
produced GRB980425. The lack of observational signatures, typically expected to be produced by 
such an off-axis jet on a 1 yr time scale, may be due to the low m < 0.1, which implies that an off-axis 
jet will become observable only on > 10 yr time scale. 

Subject headings: gamma rays: bursts and theory — supernovae: general — supernoave: individual 
(SN1998bw)— X-rays: general 



1. INTRODUCTION 

Recent Chandr a ob servations of SN1998bw 
(jKouveliotou et al~l l2004f) have allowed to identify 
one of the X-ray sources detected earlier by BeppoSAX 
(jPian ct al. 2000) as associated with the supernova. 
We discuss here the implications of the detected X-ray 
emission, in particular to the suggested association of 
SN1998bw with GRB980425. 

The association of GRBs with type Ib/c supernovae 
is motivate d by the coincidenc e of GRB980425 and 
SN1998bw fCala ma et 311 11998'). and by the identifica- 
tion of a SN1998bw-l ike spectrum in the optical after- 
glow of GRB0303292 llSta,nek et al. ll2003l Iffk^th et al. I 
|2003l). The 7-ray luminosity L^j^o. ^ lO^^erg/s of 
GRB030329, inferred from the redshift z=0.1685 of its 
host galaxy, lies within the range of t ypical cosmologi- 
cal GRBs, L^jso. ^ lO^^iigj-g/g (g.g. ISch^idt 2001) . 
The subscript "Iso." indicates luminosity derived assum- 
ing isotropic emission. The association of GRB980425 
with SN1998bw sets the distance to this burst to 38 Mpc 
(for Hq = 65km/sMpc), implying that its luminosity is 
nearly 5 orders of m agnitude lower tha n that typical for 
cosmological GRBs ijPian et al. Il200n|) . 

Two hypotheses are commonly discussed, that may ac- 
count for the orders of magnitude difference in luminos- 
ity. First, it may be that SNe Ib/c produce two different 
classes of GRBs, with two different characteristic lumi- 
nosities. It is now commonly believed that long dura- 
tion, T > 2 s, cosmological, L^jso. — lO^^erg/s, GRBs 
are produced by the collapse of SN Ib/c progenitor stars. 
It is assumed that the stellar core collapses to a black- 
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^ See, however, Katz (1994), who suggests that a SN like emis- 
sion may result from the impact of the relativistic GRB debris on 
a nearby dense gas cloud. 



hole, which accretes mass over a long period, ~ T, driv- 
ing a relativistic jet that penetrates the m antle /envelope 
and then p roduces the observed GRB (Wooslcv Il993t 
iPaczvnskil 119981 iMacFadven & Wooslev 1299) . This 
scenario is supported by the association of GRB030329 
with SN2003dh, and by additional evidence for op- 
tical sup ernovae emission in several GRB afterglows 
oom I I2OOI . The origin of a second, low-luminosity, 
class is unknown. It may be, e. g. , due to supernova 
shock break-out llColgate I Il968t iWooslev fc Weaver I 
ri98fillMatzner fc McKee 111 999HTa,n et a,l. Il200l: 

A second possibility is that GRB980425 was 
a typica l, cosmologica l GRB jet viewe d off- 
a,xis (Na,^;;;;;^ 179981 lEichler fc Leyiusoni 11995 
Wooslev. Eastman fc Schmidt I Il999t iGranot et al. I 
l2002HYamazaki. Yonetoku fc Nakamura Il2003j) . Due to 
the relativistic expansion of jet plasma , with Lorentz 
facto r F > 10 during 7-ray emission (|Krolik fc Pier I 
Il99ll IBaring I il993), 7-rays are concentrated into a 
cone of opening angle comparable to the jet opening 
angle 9j (assuming 9j > l/F). Thus, if the jet is 
viewed from a direction making an angle larger than 
9j with the jet axis, the 7-ray flux may be strongly 
suppressed. In this scenario, strong radio emission, 
~ 103° i^inr''Tj„erg/s Hz, is expec te d at ^ 1 yr delay 
fFrail . Waxman fc Kulkarni I 120001: |L ivio fc Waxman I 
2000) as the jet decelerates to sub-relativistic speed 
and its emission becon ies nearly isotro pic (for a more 
detailed discussion see iWaxmanl 120041 and references 
therein) . 

We have recently shown (Waxman 2004) that the 
low radio luminosity of SN1998bw, compared to that 
expected from a decelerated GRB jet at ~ 1 yr de- 
lay, may be consistent with the off-axis jet interpre- 
tation of GRB980425 provided that either the mag- 
netic field energy fraction behind the shock driven 
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by the jet into the wind is atypicaUy low, cb < 
10"^, or that the density of the wind surrounding 
the progenitor is lower than typically expected, rh = 
(M/10-5Afoyr-i)/(i;w/103kms"i) ~ 0.1. Lower val- 
ues of m and of es reduce the specific luminosity at the 
transition to sub-relativistic expansion. A low value of 
m <C 1 further delays the time at which the flow ap- 
proaches sub-relativistic expansion. We consider the low 
eB scenario less likely, since we expect e^, which is de- 
termined by the shock micro-physics, to be similar for 
different bursts, and es close to equipartition is inferred 
from radio observations of other bursts. The latter sce- 
nario, rh ~ 0.1, is consistent with the constraints im- 
posed on rh by the observed ra dio emission from the su- 
pernova (Kulk arni et al. II199 8'). which is interpreted as 
due to synchrotron emission from a shock driven into the 
wind by a mild ly relativistic shell ejected by the super- 
nova explosion iKulkarni et al. |ll998t IWaxman fc Loebl 
[1999; Chevalier fcLilll999D . In this scenario, transition 
to sub-relativistic expansion is expected over ~ 10 yr 
time scale (, Waxman . 20041 . 

The recent identification of the X-ray emission of 
SN1998bw was argued to support the interpretation that 
GRB980425 was an off- axis "standard" GRB jet asso- 
ciated with SN1998bw ijKouveliotou et al. ll200^ . We 
show here that the observed X-ray emission may be nat- 
urally explained as synchrotron emission from a shock 
driven into the wind surrounding SN1998bw progenitor 
by the mildly relativistic shell that was inferred to ex- 
ist from earlier radio observations. In 321 we construct a 
simple model that describes the emission from a shock 
driven into the wind by a s ub-rel ativistic shell, based on 
the analysis of IWaxmanI l)2004j) . In we show that 
the model can account for both X-ray and radio observa- 
tions. While X-ray observations, or radio observations, 
alone do not allow to determine all model parameters, 
combined X-ray and radio data over constrain the model, 
and therefore enhance our confidence in its validity. Our 
conclusions are summarized in 21 

2. MODEL 

2.1. Dynamics 

Let us consider a sub-relativistic shell ejected by the 
supernova explosion, with mass M, total energy E and 
initial velocity v — f3c. We assume the shell's thermal 
energy is dominated by its kinetic energy, since the ther- 
mal energy deposited in the shell by the passage of the 
supernova shock leads to expansion and post-shock ac- 
celeration which converts the thermal energy to ki netic 
energy (As pointed out bv lMatzner fc McKee 1119991 this 
post-shock acceleration is essential for achieving mildly 
relativistic shell velocity). The shell is assumed to prop- 
agate into a density profile created by stellar mass 
loss. We assume a density profile 
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where M is the mass loss rate and Ww is the wind velocity. 

As the shell expands, it drives a shock wave into the 
surrounding wind, and a reverse shock is driven back- 
ward (in the shell frame) into the shell. It is straight- 
forward to show that as long as the mass of shocked 
wind plasma, AIw{r) — AnKr, is small compared to the 



shell's mass, the shock driven backward into the shell 
does not lead to significant deceleration. Denoting by 
u < V the velocity of shocked wind (and shell) plasma, 
pressure balance between the forward and reverse shocks 
implies PwU^ = Ps(y — w)^, where pw = M^/Airr'^ and 
Ps = rjsM I Airr^ denote wind and shell density respec- 
tively. Here r/rjs is the thickness of the ejected shell. 
For < M we find {v -u)/v = {M^/r]sM)^^'^ . Thus, 
at small radii, where Af^ ^ M, the shell expands with 
a time independent velocity, v = (3c. At this stage, the 
velocity of the shock propagating into the wind is also 
time independent and given by (7 + l)w/2, where 7 is 
the adiabatic index of the wind plasma. 

Significant deceleration of the shell begins at a radius 
where the shocked wind mass, M^(r), is similar to the 
shell's mass M . At larger radii the flow approaches the 
Sedov-von Neumann- Taylor solutions describing expan- 
sion of a spherical strong shock wave into a density 
profile (e.g. Chapter XII of Zel'dovich & Raizer 2002). 
In these solutions, the shock radius is given by 



1/3 



(2) 



^(7) is a dimensionless parameter of order unity, which 
may be approximated by ^(7) — (3/2) [(7 -I- 1)^(7 — 
l)/27r(77 - 5)]i/3 llWaxmanll200l . For 7 = 5/3, ap- 
propriate for sub-relativistic flow, we have ^ = 0.73. At 
the radius where Myj{r) ~ M, the velocity given by the 
Sedov-von Neumann- Taylor for the appropriate energy E 
is similar to the initial shell velocity v. We therefore de- 
fine the deceleration radius, i?dcci as the radius at which 
the Sedov-von Neumann- Taylor post-shock fiuid velocity 
equals v, 2R/{'^ +1) = v, 
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M/v„ 



2.1 X 10 



15 -^49 



(3'^rh 



(3) 

where E — 10'*^£'49 erg. The deceleration time, idee, is 
defined as 



Rd 



(7 + l)v/2 



0.62|^day. 



(4) 



For t <^ tdcc we therefore expect shell expansion at 
constant velocity, v, and a shock driven into the wind at 
a constant speed, (7 -I- l)v/2. For t ^ ^dec, we expect 
the expansion to follow the Sedov-von Neumann- Taylor 
solution, Eq. |21l, with shock velocity R oc t~^/^. 

2.2. Synchrotron emission 

Let us consider synchrotron emission from electrons 
accelerated to relativistic energies by the collision- 
less shock driven into the wind. We assume that 
a constant fraction Cg (cs) of the post-shock ther- 
mal energy is carried by relativistic electrons (mag- 
netic field), and that the electron distribution func- 
tion follows a power-law, dlnne/dln7e = —p for 
7e > 7m- In what follows, we assume p = 2, 
as ob served for non-relativistic l)Blandford fc EichlerJ 
1987 1 as well as for relativistic shocks (Waxma^ 



1997J). This power law is prod uced by Fermi accelera- 
tion in coUisionless shocks (iBlandford fc Eichler I [1987t 



tion in coUisionless sh ocks (iBlandford fc Eichler J1983 
iBednarz fc Ostrowski I 119981: lAchterberg et al. I l2001(l 
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although a first principles understanding of the process 
is not yet available. 

At times t 3> tdcc the flow is well described by the 
spherical non-relativistic self-similar solution, where the 
shock radius is given by Eq. For simplicity, we as- 
sume that at this stage the shocked plasma is concen- 
trated into a thin shell behind the shock, R/rj ^ R, 
within which the plasma conditions are uniform (the 
Chcrnyi-Kompaneets approximation). This implies, in 
particular, ?7 = (7 + l)/(7 — !)■ The synchrotron emis- 
sion obtained under the above assumptions was derived 
in ijWaxman 112004^) . The scaling of magneti c field am- 
plitud e B and of 7™ with time, at i » tdcc, is (jWaxman I 

l2?inl 

Boct-\ 7„oct-2/3^ (5) 
and for p = 2 the specific luminosity is 



-3/2 



VlkeV; VlOOday; 

Eq. ^ holds at frequencies well above the characteris- 
tic synchrotron frequency of the lowest energy electrons 
(7e — 7m), and below the cooling frequency, Vc- Vc, 
the characteristic synchrotron frequency of electrons for 
which the synchrotron cooling time is comparable to the 
shock expansion time, t, is 



2 X 10'''(3eB) 



/ t 



V lOOday 



Hz. 



(7) 



Let us now consider synchrotron emission from the 
shock driven into the wind at times t <C tdcc- At frequen- 
cies v > i^c, the specific luminosity (for p = 2) is given by 
= Lm{'^c/'^m)~^^'^{'^/'^c)~^ = Lm{vcVrnY/'^ /v. Here, 



is the specific intensity at v 



the character- 



istic synchrotron frequency of the lowest energy elec- 
trons (7e = 7m)- In order to determine the time de- 
pendence of Lj/, we first note that Lm is proportional 
to the product of the number of radiating electrons, Nf., 
and the magnetic field strength, B, Lm oc NeB, and that 
Vm oc 7^-0- The Lorentz factor 7c of electrons that cool 
on a time ~ t is given by comparing the synchrotron cool- 
ing time, tsyn oc Xj^Ji'^ to t, yielding 7^ oc 1/tB^ and 

oc 7^5 oc t^^B~^. The dependence of i^rm '^cj and 
on time is now obtained by noting that the time inde- 
pendent shock velocity implies time independent thermal 
energy per shocked particle, i.e. 7™ oc t'^, thermal energy 
density oc oc (since the density drops as l/r^), i.e. 
B oc t"^, and oc t. Thus, L,„ oc t'^, Vm oc 7^-6 t~^, 
Vc <yit and Ly oc t'^ for v > Vc- Synchrotron emission at 
times t tdcc is therefore described by 



(8) 



with given by Eq. Q. Combined with Eqs. Q 
and ©, this determines the synchrotron emission pre- 
dicted by the model at t < idcc- 

The following line of arguments shows that synchrotron 
emission from the shock driven backward into the ex- 
panding shell, which exists at t < tdcc, is dominated by 
the emission from the shock driven into the wind, under 
the assumption that and are similar for the post- 
shock plasma behind the two shocks. Since the thermal 



energy density is the same in the shocked wind and shell 
plasmas (pressure equilibrium) , the magnetic field is the 
same in the two regions (assuming similar es for both 
shocks) and the product of the electron number den- 
sity and their minimum Lorentz factor, ricjrin which is 
proportional to the thermal energy density, is the same 
in the shocked wind and shell plasmas. Similar B im- 
plies similar for both shocks, and hence the ratio of 
the specific synchrotron luminosity of the shocked shell 
plasma and shocked wind plasma, L^^s/L^, is L^^g/L^, = 

Lm.sVm^s/LmVm'^ = Ne,sjm,s / Nejm = Ag/A, whcrC the 

subscript s denotes quantities related to the shocked 
shell plasma, and A (A^) is the width of shock wind 
(shell) plasma. The thickness of the region of shocked 
wind/shell plasma is proportional to the velocity of the 
shock propagating into the wind/shell plasma, measured 
in the shocked plasma rest frame, A^/A = (v — u)/u [see 
discussion following Eq. (^]. Thus, for ^ M we find 
L,.s/L, = {M^/r]sM)^/^ « 1. 

The following point should be emphasized here. The 
model described above is highly simplified. For exam- 
ple, the wind density is assumed to follow a pure l/r^ 
law and to be free of inhomogeneities, and the expanding 
shell is assumed uniform with no internal (density, veloc- 
ity) gradients. In comparing model predictions with ob- 
servations one should therefore expect an approximate, 
rather than exact, agreement. Construction of a more de- 
tailed, and more accurate, model requires a more detailed 
knowledge of, for example, the wind density distribution 
and its homogeneities and of the internal structure of 
the shell. We find that such a more detailed analysis is 
not warranted by the present data, since the number of 
free parameters in a more detailed model would be larger 
than the observational constraints. 

3. IMPLICATIONS TO SN1998BW 

The luminosity of the X- ray source associated wit h 
SN1998bw, given in Fig. 3 of lKouveliotou et al. I l|2004j) . 

may be approximately described as 



L 



X 



6 X 10 erg/s x 



1, t <110 day; 

(t/110day)"3/2, otherwise. 

(9) 

The t 3/2 decline at late times is consistent with the 
prediction of Eq. ©, describing synchrotron emission 
from the shock driven into the wind at the self-similar 
stage of expansion, t > tdcc [see Eq. Q], at frequen- 
cies lower than the cooling frequency, given by Eq. 0. 
The time independent flux at earlier time is consistent 
with synchrotron emission from the shock driven into 
the wind at the stage of expansion at constant speed, 
t < tdoc, at frequencies higher than the cooling fre- 
quency [see Eq. JSJ]. The observed X-ray luminosity 
may therefore be interpreted as due to the shock driven 
by a sub-relativistic shell into the wind, for which de- 
celeration starts at tdcc ~ HO day and the cooling fre- 
quency passes through the X-ray band at t ~ tdoc, i-e- 
/iiyc(t = llOday) « 1 keV. 

Using Eq. (0, the requirement hi/dt — 110 day) ~ 
1 keV implies 



(3eBm)3/4 « 10- 



(10) 



Using this result in Eq. ((BJ, the observed luminosity at 



4 



Waxman 



t = tdoc ~ 110 implies 



6. 



(11) 



Using this result in Eq. Q), the requirement tdec ~ 
110 day implies 

/J^m « 0.03/(3ee). (12) 

Our model for the X-ray emission of SN1998bw is 
consistent with its radio emission, which is also inter- 
preted as synchrotron emission from a shock wave driven 
into the wind by a mildly relativistic shell ejected by 
the supernova explosion. Within the framework of this 
model, radio observations imply electron energy frac- 
tion close to equipartit ion, See ^ 1 (see Eq. (17) of 
IWaxman &: Loeb1ll999(l . Near equipartition value of 
is also inferred from radio observations of GRB970508 for 
the shock driven by the fireball into surrounding plasma 
at late stages, where the shock becomes mildly rela- 
tivistic l,Frail. W axman & Kulkarni 2000) . We therefore 
adopt Bee ~ 1, which implies, using Eq. Hll(l . £^49 ss 6. 

Radio observations further allow to determine the val- 
ues of es, /3 and E as functions of m. The lowest al- 
lowed value of m (and of E) is that corresponding to 
magnetic field equipartition, with larger values of m 
implying lower values of cb (an d larger values of E). 
Model predictions were shown dKulkar ni et al. 1 119981: 
IWaxman fc Loeb ll999trChevalier fc Li I1999D to be con- 
sistent with radio observations for, e.g., {m = 0.04, 
es « 0.1, Eig w 3, /3 w 0.8}, and for {m = 6, Eig « 50, 
cb ~ 10^^, P w 0.6}. The energy inferred from X-ray 
observations, Eq. (|11|) . implies that the former choice of 
parameters, {m = 0.04, eb ~ 0.1, -E49 w 3} with (3 « 0.8, 
should be chosen for consistency with observations. 

X-ray and radio observations therefore over constrain 
the model: All model parameters are determined from 
radio observations with the additional X-ray constraint 
given by Eq. (|ll|l . The inferred parameter values are con- 
sistent with the two additional independent X-ray con- 
straints, Eqs. (O and l(T^ . 



4. DISCUSSION 

Simple analytic expressions are given in § |21 for the 
specific luminosity emitted by a shock wave driven into 
a wind surrounding a supernova progenitor by a sub- 
relativistic shell ejected by the supernova explosion. We 
have shown in §|21that the observed X-ray and radio emis- 
sion from SN1998bw are consistent with such a model. 
In particular, the approximately time independent X-ray 
fiux at early time, t < 100 day, and the t~^/^ decline 
of the flux at late time [see Eq. are characteristic 
for the emission from a shock driven into a wind by a 
shell which suffers significant deceleration at t > 100 day 
[see Eqs. ©, ©]. Combined X-ray and radio data al- 
low to determine all model parameters: The shell ki- 
netic energy is £' « lO'*^'^ erg and its initial velocity 
is /3 « 0.8; The normalized mass loss rate of the wind, 
m = (M/10-^MQyr-i)/(ww/10^knis-i), is m « 0.04; 
The post-shock magnetic field is not far below equiparti- 
tion, eB ~ 0.1. Combined X-ray and radio data over con- 
strain the model, and therefore enhance our confidence in 
its validity. We note, that the inferred energy in mildly 
relativistic ejecta, T(3 > 1 where F = (1 — /3^)~^/^, is an 
order of magnitude higher than the value obtained from 
analyses of shock emergence fr om the assumed p rogeni- 
tor of SN1998bw, ~ lO'*^-^ erg (iTan et al. 1127)01 . 

X-ray observations do not provide therefore evidence 
for existence of an off-axis "standard" GRB jet associated 
with SN1998bw, that may have produced GRB980425. 
The lack of observational signatures typically expected 
to be produced by such an off-axis jet on a 1 yr time scale 
may be due, however, to the low value of m, m < 0.1, 
which implies that an off-a xis jet will becom e observable 
only on > 10 yr time scale l|Waxman Il2004j) . 
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